Coexistence of color superconductivity and 
chiral symmetry breaking within the NJL 

model 

D.Blaschke 1 * 2 , M.K.Volkov 2 , and V.L.Yudichev 2 

Fachbereich Physik, Universitat Rostock, D-18051 Rostock, Germany 
2 ' Joint Institute for Nuclear Research, 141980 Dubna, Russian Federation 

Abstract 

The phase diagram for quark matter is investigated within a simple 
Nambu-Jona-Lasinio model without vector correlations. It is found 
that the phase structure in the temperature-density plane depends 
sensitively on the parametrization of the model. We present two 
schemes of parametrization of the model where within the first one 
a first order phase transition from a phase with broken chiral symme- 
try to a color superconducting phase for temperatures below the triple 
point at Tt = 55 MeV occurs whereas for the second one a second or- 
der phase transition for temperatures below Tt = 7 MeV is found. 
In the latter case, there is also a coexistence phase of broken chiral 
symmetry with color superconductivity, which is a new finding within 
this class of models. Possible consequences for the phenomenology of 
the QCD phase transition at high baryon densities are discussed. 
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1 Introduction 



The phenomenon of color superconductivity [1112113111130111013 UH1 

nu na nsi im chi nni czi na cni 1201 1211 iza 123 1211 12a iini is of general 

interest, in particular, in studies of the QCD phase structure [01 EH EH E3 
IT71 HH1 EE Em EH E21 EH] and applications in the astrophysics of compact 
stars [IH1EZI- Observable consequences are expected for, e. g., the cooling 
behavior [2E1E0I- Different aspects have been investigated so far, whereby 
models of the NJL type have been widely employed [301 EH E21 EH] in studies 
of the phase structure in the vicinity of the hadronization transition. 

Recently, it has been shown in these investigations that for low temper- 
atures (T) and not too large chemical potentials (//) the two-Flavor Color 
Superconductivity (2SC) phase is favored over alternative color supercon- 
ducting phases [201 EH EH ES] • According to [213, the Color-Flavor-Locked 
(CFL) phase occurs only at \i > 430 MeV. 

It is generally agreed that at low temperatures the transition of the mat- 
ter from the phase with broken chiral symmetry to the color superconducting 
phase is of the first order (see e. g. P2EJ). From the point of view of phe- 
nomenological applications, as e.g. in compact star physics, the order of 
the phase transition to quark superconducting matter plays an important 
role. The conclusion about the first order phase transition was drawn within 
models without vector interaction channels taken into account; the vector 
interaction has been considered in few papers [HH E2] • It was found that the 
presence of quark interaction in the vector channel moves the critical line in 
the in /i — T plain to larger // [TH EH1 EI] • Recently it has been demonstrated 
[3*5] that the critical line of first order phase transition in the /i — T plane can 
have a second end-point at low temperatures, besides the well known one at 
high temperatures. The latter one could even be subject to experimental ver- 
ification in heavy-ion collisions [H] whereas the former could be of relevance 
for neutron stars. While in Ref. [HH] this feature of the phase diagram was 
a consequence of the presence of interaction in the vector channel, we would 
like to investigate in the present work the sensitivity of the phase diagram 
to the choice of model parameters without interaction in the vector channel. 
We will demonstrate that in the absence of the vector channel interaction the 
phase transition is not necessarily of the first order, thus revising statements 
in Refs. |33 EH - 

It is worth noting that some progress has recently been done in lattice 
calculations. There are methods being developed that allow to extend lattice 
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results to the case of finite chemical potentials |4TH 141[ However, these 
methods are valid only for small chemical potentials (see e. g. |40j), below the 
conditions at which the color superconductivity phase is expected to form. 

The structure of our paper is as follows. In Sect. 2, a chiral quark model 
is introduced, its Lagrangian is given and the model parameters are fixed 
from the vacuum state in two different schemes. Temperature and chemical 
potential are introduced into the quark model in Sect. 3, using the Matsubara 
formalism. The conclusions and a discussion of the obtained results are given 
in Sect. 4. 

2 NJL model with the scalar diquark channel 

In order to study the quark matter phase diagram including color supercon- 
ductivity, one should generalize the concept of the single order parameter 
related to the quark-antiquark condensate in the case of chiral symmetry 
breaking to a set of order parameters when condensation can occur in other 
interaction channels too. The simplest extension is the scalar diquark con- 
densate (i/jip) for u and d quarks 



which is an order parameter characterizing the domain where the color sym- 
metry is spontaneously broken and the quark matter finds itself in the (two- 
flavor) color superconducting (2SC) state. This quantity is the most im- 
portant one among other possible condensates that can be constructed in 
accordance with the Pauli principle j2S]. In the matrix C is the charge 
conjugation matrix operator for fermions 



The matrices t 2 and A2 are Pauli and Gell-Mann matrices, respectively. The 
first one acts on the flavor indices of spinors while the second one acts in the 
color space. 

If the electroweak interaction is discarded and only the strong coupling is 
in focus, the resulting quark matter phase diagram is essentially determined 
by nonperturbative features of the QCD vacuum state. One therefore has 
to resort to nonperturbative approaches to describe the behavior of parti- 
cles at various conditions, ranging from cold and dilute matter up to the 



5 = {ipi^r 2 X2Cii T ) 



(1) 



C = i7o7 2 - 



(2) 
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hot and dense one. A reliable and widely tested model to nonperturbative 
strong coupling QCD is provided by the Dyson-Schwinger equations 
however, for qualitative studies like the one we attempt here it proves to 
be too complex. Therefore, we will use here a simple and tractable nonper- 
turbative model of quark interaction, the Nambu-Jona-Lasinio (NJL) model 
[3TH EU E21 CUE EH ESI j which has been extensively exploited for the descrip- 
tion of the properties of the light meson sector of QCD (also to describe 
the color superconductivity phase EH1 EZI) and proved to be a model 
respecting the low-energy theorems. Before we proceed to the case of finite 
temperature and density, the model parameters that determine the quark 
interaction should be fixed. This shall be done for the vacuum state where 
hadronic properties are known. We will assume, according to common wis- 
dom that, once fixed, these parameters (originating from the nonperturbative 
gluon sector of QCD) will not change, even in the vicinity of the transition 
to the quark matter. This transition is thus caused by medium effects in the 
quark sector only. 

2.1 Lagrangian 

In the present paper we restrict ourselves to the two-flavor case, leaving 
the strange quark and effects related to it beyond our consideration. As 
we constrain ourselves to only two order parameters, the quark and scalar 
diquark condensates, during our investigation, the interaction of quarks will 
be represented in the Lagrangian by SU(2)l x SU(2)r symmetric scalar, 
pseudoscalar quark- ant iquark, and scalar diquark vertices: 

C = ^ - m°)tp + C qi! + C qq , (3) 

As = % KW 2 + (WV) 2 ] , (4) 

£qq = ^{^1^2^2C^ T )(i) T Ci^T2\2^), (5) 

where m° is the diagonal current quark mass matrix m° = diag(m u , ra^), 
G and H are constants describing the interaction of quarks in the scalar, 
pseudoscalar, and scalar diquark channels, respectively. We work in the 
isospin symmetric case (mj = m? d = m°), thus m° = m°lf. 

By the standard Hubbard- Stratonovich procedure , we introduce auxil- 
iary scalar (ex), pseudoscalar triplet (7?), and diquark (A, A*) fields together 
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with Yukawa-like terms in the Lagrangian density instead of the four-quark 
vertices: 



C = $ — m° + cr + ij^TTr)^ — 



2G 2H 

In order to integrate out the quark degrees of freedom by Gaussian path 
integration, it is appropriate to represent the quark fields by the bispinor 



and to introduce the matrix propagator S(p): 

Integrating over q(x) and we then obtain an effective Lagrangian in 

terms of collective scalar and pseudoscalar quark-antiquark and scalar di- 
quark excitations. Here, we restrict ourselves to the mean-field approxi- 
mation, leaving the next-to-leading order corrections in the 1/N C expansion 
beyond our model. Finally, the effective Lagrangian density reads 

a 2 + 7f 2 |A| 2 r d^p 1 

The trace in (jHJ) is taken in the Dirac, color, and flavor space. The matrix 
M contains a and n fields: 

M = (m° - a)l - i l5 r a Tr a , (10) 

the sum over a = 1, 2, 3 is assumed, and 1 = l c • 1/ • Id. 

As it was mentioned above, we are working in the mean-field approxima- 
tion and the quark condensates are of interest. Therefore, further study can 
be performed in terms of the effective potential 

Vcs = - lim — / d 4 xC cS (11) 
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where V4 is 4-dimensional volume. The vacuum expectation values of the 
collective variables a, n, A, and A* determine the absolute minimum of V e g. 
They are given by the equation 

dVcB = dV cS = dVeg _ 

da ~ OA ~ 8A* ~ ' 1 ' 

A priori it is known that in the vacuum only the a field acquires a nonva- 
nishing expectation value. The diquark fields A, A* are expected to have 
nonzero mean values only in dense matter. The mean value of the pseu- 
doscalar isotriplet field n is always equal to zero, therefore we omit it here- 
after. 

Having solved Eq. (fT2j) for the field a, one can work in terms of the 
constituent quark mass m, connected with the current quark mass by the 
gap equation 

m -m= (a) = 2G(ip^) . (13) 

In the chiral limit (m° = 0), the constituent quark mass is proportional to 
the quark condensate and thus can be treated as the order parameter. 
In the NJL model the quark condensate is 

= -4mlf(m), (14) 

where 

xw-T&fa'-ftirh- (15) 

The divergence in I^{m) is eliminated by means of a sharp 3D cut-off at the 
scale A. 



2.2 Parameter fixing 

In our model we have four parameters: the four-quark interaction constants 
G and H, cut-off A, and the current quark mass m°. Without diquarks, there 
are only three: G, A, and m°. They are fixed by the following relations: 

1. The Goldberger-Treiman relation (GTR): 

m = g n F n , (16) 

where _F^ xp pa 93 MeV is the pion weak coupling constant and g n de- 
scribes the coupling of a pion with quarks g n 7iipTip 

-2 ,,a,_, ,a,^> ~' N c fe(A*-p)d<p 
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2-a. The quark condensate (QC) from QCD sum rules 



(^V')qcdsr 



4mlf(m) « (-240 MeV) 3 . 



(18) 



2-b. The decay constant g p for the p — > 2tt (R2PD) process 



g„ = y/6g v , g? p « 6.1 . 



(19) 



The 7r — ai transitions are omitted here. 



3. The current quark mass m° is fixed from the GMOR relation: 



Ml 

7T 



2m° (ipip) 



M exp 

7T 



140 MeV . 



(20) 



F 2 

7T 



In the chiral limit M n = 0, m = 0. 

4. With the diquark channel included, there is an additional parameter 
H which can be fixed as H — 3/4G from the Fierz transformation (as 



In the item 2, we have given two alternatives: one can either use the 
value of the quark condensate taken from QCD sum rule estimates or demand 
from the model that it should describe the p — > 2ir decay. The latter is well 
observable in experiment contrary to the quark condensate. 

For simplicity, we perform all calculations in the chiral limit m° = 0. 
In this case, when investigating the hot and dense quark matter, the bor- 
ders between phases turn out to be sharp and the critical temperature and 
chemical potential are well defined. With the finite current quark mass, the 
transitions from one phase to the other become smooth. 

As a result, one obtains two different parameter sets shown in Tabled 

In the Type I parameter set the interaction of quarks is stronger, the UV 
cut-off is smaller, and the constituent quark mass is greater. One can calcu- 
late the dimensionless constant GA 2 . It equals 4.6 for the Type I and 3.72 for 
the Type II, respectively. As we will see further, these two parametrizations 
result in qualitatively different phase diagrams. 

1 Somc authors use H = 1/2G. It turned out that within our model the resulting phase 
diagram is not much affected if one makes the choice in favor of H = 1/2G. However, it 
would be preferable to fix the constant H from some observable, e. g. from the nucleon 
mass. 



e. g. inj^Ol) \ 
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3 NJL model at finite T and fi 
3.1 Thermodynamical potential 

We extend the NJL model to the case of finite temperatures T and chemical 
potentials /x, applying the Matsubara formalism, and restrict ourselves to 
the isospin symmetric case where up and down quark chemical potentials 
coincide. The thermodynamical potential per volume is 



fi(3» 



t— -Trln ( — S 1 (iu„,p) 
(2tt) 3 2 \T K ,FJ 



+ 



IAI 



2G • 2H ' (21) 
where u; n = (2n + l)ixT are Matsubara frequencies for fermions, and the 
chemical potential is included into the definition of inverse quark propagator 



S (po,p) 



i> - M - /J70 A75T2A2 
-A*7 5 r 2 A 2 i> - M + /i7o 
The expression in (12 lj) can be simplified using the equations 

II III ! . 1 I /.J.', ./ill = I I t I I ■ 



(22) 



+ 2 



-2w y , _2, 



In 



^4 



(23) 



and 



n=— 00 

= - 



2^2 



£ ± + 2Tln[l + exp(-£ ± /T)]. (24) 



1 + exp 



T 



+ 2Tln 



1 + exp ( - — 



+ 4 ( (£+ + 6(A 2 -f) +2Tln 



+ 2Tln 



1 + exp 



T 



1 + exp 

m 2 |A| 2 
+ 2G + 2F 



T 



(25) 



where 

e = a/jo 2 + m 2 , e ± = e ± fi, 



(26) 



£± = ^/PF+W- (27) 
The cold matter limit T = looks as follows: 

n(o,*0 = -/^[2(| e + | + | e -|) + 4(E + + J E-)] 

m 2 IAI 2 

The thermodynamical potential cannot be calculated in closed form for ar- 
bitrary T and fi. However, in the cold matter limit one can easily obtain 
analytic expressions for the thermodynamical potential or its derivatives if 
only one of the collective variables o or |A| has a nonvanishing average value. 
This allows to find what kind of phase transition is to be expected for different 
parameter choices. 

We evaluate the remaining 3D momentum integrals numerically and cal- 
culate the value of thermodynamical potential at different T and fi for the 
two types of model parameter sets. The equilibrium state for each T and 
fi is determined by (a) = —in and (|A|) corresponding to the minimum of 
<>('/•.//). 

3.2 Numerical results: Type I 

It is quite illustrative to look at the contour plots of the thermodynamical 
potential. For several values of fi at T = they are shown in Figs. where 
one can follow the appearance and disappearance of local minima, maxima, 
and saddle points of the thermodynamical potential with increasing chemical 
potential. For zero temperature and chemical potential we have, as expected, 
a nonzero constituent quark mass (quark condensate) corresponding to the 
absolute minimum of the thermodynamical potential at m ~ 350 MeV and 
(|A|) = in Fig. [I] At a certain chemical potential, a new local minimum 
related to the diquark condensate near (|A|) ~1 10 MeV and m = (Fig. EJ), 
but it does not yet give the absolute minimum. There is also a local max- 
imum around m ~ 200 MeV and (|A|) = 0. As the matter becomes more 
dense, the second minimum lowers until it becomes degenerate with the first 
minimum while the average value of a (or —m) remains almost unchanged 
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(see Fig. |2J). Above the corresponding (critical) chemical potential /i c ~ 321 
MeV, the second minimum becomes the absolute one and a first order tran- 
sition occurs, during which (er) discontinuously changes to zero while the 
diquark condensate acquires nonzero value breaking the color symmetry of 
the strong interaction. This characterizes the color superconducting phase 
transition in quark matter. Furthermore, the local minimum on the m axis 
merges the saddle point (see Fig. |3J) and, at still higher /x, only the local 
minimum on the A axis near |A| ~ 130 MeV and m = remains. 

At a fixed chemical potential above /i c , with the temperature rising, the 
average value of |A| decreases until it reaches zero at the critical temperature 
T c which can be roughly estimated using the BCS theory formula 



Above this temperature quark matter is in the symmetric phase 2 where 
the chiral and color symmetries are restored. Finally, we obtain the phase 
diagram shown on Fig. El with three phases: the hadron phase, 2SC phase, 
and symmetric phase. All three phases coexist at the triple point: T t ~ 55 
MeV and fj, t « 305 MeV. 

3.3 Numerical results: Type II 

As for the Type I parameter set, at zero T and fi only the constituent quark 
mass m, being the order parameter for the chiral condensate, is nonzero, 
whereas the diquark gap A vanishes. However, the vacuum value of m is 
lower than that for the Type I and, with the chemical potential increasing, \i 
becomes equal to the vacuum value of m before the second local minimum, 
corresponding to the diquark condensate, appears. At further increase of ji 
the constituent quark mass decreases, and it would vanish at \i = fix, 



if the diquark condensate did not appear. Actually, at the critical value \L\ 
both the quark condensate and the diquark condensate are small but nonzero. 

The changes of the local extrema for increasing chemical potential are 
similar to those shown in Figs. EHU The cases of dilute (/i = 0) and very 

2 According to recent investigations |48j . a so-called pseudo-gap phase as a precursor of 
color superconductivity can occur in this region. 




(29) 




(30) 
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dense matter (// = 400 MeV) are qualitatively analogous, only the absolute 
values of m and |A| at which the local minima are found are different. At 
intermediate densities, however, there is a qualitative difference. Within a 
very narrow range of values of the chemical potential, there exists a new 
phase of massive superconducting matter. One can see this in Fig. El for 
fi = 286 MeV. At higher [i the chiral symmetry is restored and the quark 
matter is in the pure superconducting phase. A possibility of the chiral 
diquark condensates to coexist at certain condition has been already noticed 
in Ref. [H] 

Thus, for the Type II parameter set, the transition from the hadronic to 
the superconducting phase is of the second order. In this case there are no 
degenerate local minima in the thermodynamical potential separated by a 
barrier. This behavior is unlike to what is commonly expected for a cold and 
dense matter but it parallels the findings of Ref. jHH] where vector interactions 
are responsible for this behavior. 

The average value of |A| is much smaller than for the Type I parameter 
set. As a consequence, the border between 2SC and the symmetric phases 
of quark matter lies at noticeably lower temperatures. The phase diagram 
obtained in our model for the Type II parameter set is shown in Fig. [7| 

4 Conclusion 

In the framework of the simple NJL model for two flavors, a phase diagram 
is obtained for T = 0-200 MeV and fi = 0-450 MeV. Three phases are found 
for the Type I parameter set and four phases for the Type II parameter 
set. The critical temperature and chemical potential obtained in the Type I 
scheme differ from those obtained with the Type II parameter set. At T = 0, 
/i c w 320 MeV for the Type I parameter set and fi c ~ 288 MeV for the 
Type II. The corresponding quark densities differ by a factor 1.5 - 1.7. The 
critical temperature for the Type II parameter set is as low as 7 MeV and 
thus much closer to critical temperatures for the paring instability in nuclear 
matter systems (see j5U]) whereas for the Type I parameter set the critical 
temperatures are an order of magnitude larger. This striking difference in the 
critical parameters obtained within the same model calls for a more detailed 
investigation of the question of model parametrization. 

In our work, the constant H was not obtained from a fit to observable 
data. Instead, Fierz transformation arguments have been used to fix the 
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ratio H/G = 3/4. A parameterization would be favourable where (in the 
spirit of the Type II model) experimentally measured quantities, like the p 
meson width, are used rather than non-observable ones (quark condensate 
etc.). It would therefore be more consistent to fit the constant H from 
baryon properties, see [HH E2] and also to go beyond the mean field level of 
description. These investigations shall be performed in future work where it 
remains to be clarified which critical parameters for the color superconducting 
phase transition can be considered more realistic and of which order the phase 
transition is. 
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Tables 





A 
[GeV] 


G 

[GeV- 2 ] 


H 

[GeV- 2 ] 


m 
[MeV] 


^- (W>) 
[MeV] 


9 P 


GTR+QC 
[Type I] 


600 


12.8 


9.6 


350 


240 


9.2 


GTR+R2PD 
[Type II] 


856 


5.1 


3.8 


233 


284 


6.1 



Table 1: The model parameters for two different schemes of parameter fixing. 
The first row corresponds to the case where the quark condensate value is 
used, while in the second row the parameters correspond to the case where 
the decay width of the process p — > 2n is used. 
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Figures 





Figure 2: The contour plot for the thermodynamic potential as a function of 
of m and A at zero temperature and the chemical potential \i = 340 MeV. 
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Figure 3: The contour plot for the thermodynamic potential as a function of 
of m and A at zero temperature and the chemical potential \i =350 MeV. 
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Figure 4: The contour plot for the thermodynamic potential as a function of 
of m and A at zero temperature and the chemical potential \i = 400 MeV. 
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Figure 5: The quark matter phase diagram from the NJL model with the 
Type I parameter set. In phase I both chiral and diquark condensates vanish; 
in phase II the chiral symmetry is broken; in phase III the chiral symmetry 
is restored while the diquark condensate is nonzero. 
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Figure 6: The contour plot for the thermodynamic potential as a function of 
of m and A at zero temperature and the chemical potential \x = 286 MeV 
for the Type II parameter set. 
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Figure 7: The quark matter phase diagram from the NJL model with the 
Type II parameter set. In phase I both chiral and di quark condensates vanish; 
in phase II the chiral symmetry is broken; in phase III the chiral symmetry 
is restored while the diquark condensate is nonzero; in phase IV both the 
chiral and diquark condensates coexist. 
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